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SUMMARY 


The nonlinear beam equation is integrated numerically in 
a direct -fashion to obtain results for large deflections of 
cantilevers under tip loads of arbitrary direction. A short 
BASIC computer program for performing this integration is 
presented. Results for selected load cases are presented. 

The numerical process is performed rapidly on a modern > 

microcomputer , and comparisons in the paper with results from 1 

closed form solutions from the literature show that the 
process is accurate. j 

1 

/ 

INTRODUCTION j 

f 

j 

The engineering theory of deflection of beams is linear ) 

if. deflections are small with respect to beam length. j 

Solutions to linear beam equations are presented in 
textbooks. On the other hand, if the loading and stiffness 
of the beam are related in such a way as to lead to 
deflections not small with respect to beam length; the beam 
equations are nonlinear, and solutions are more difficult. 

The purpose of this paper is to consider the large deflection 
of a cantilever beam with an arbitrarily directed tip load. 

This particular type of beam and loading could be an element 
of a larger structural or mechanical system, and appropriate 
solutions could facilitate treatment of the beam element as a 
nonlinear spring in the- context of the large system. It 
turns out that exact, closed-form solutions to this problem 
exist in the literature (ref. 1). These solutions, however, 
are in terms of elliptic functions, and sometimes the desired 
unknowns are arguments of the elliptic functions. To obtain 
numerical results, therefore, may require iteration or 
trial-and-error calculations using tables of elliptic 
functions. In this paper, a different approach is taken. 

The nonlinear beam equation is integrated numerically in a 
direct fashion. The equations for performing this numerical 
integration are presented, and results are presented for 
selected cases. Finally, a short BASIC-1 anguage computer 
program is presented which performs the numerical integration 
of the beam equations. The program is written in Applesoft 
BASIC for an Apple 11+ microcomputer with 64K memory. 


1 

4=~ 


ht<7*T- 



SYMBOLS 


A 


El 

F 

h 

L 

M 


P 


Q 


LJ 

V 


X , Y 

XI , Y1 
XL, YL 

0 


=Fl?/EI, Dimensionless 1 oad-st i f -f ness parameter 
Beam bending stiffness 
Magnitude of tip load 

Increment in dimensionless length parameter, s 
Length of beam 
Bending moment in beam 
Load along neutral axis of beam 
Transverse shear load in beam 
=dv^/ds 
= 

Coordinate along neutral axis of beam 
= s/L 

Beam coordinates with origin at root 
Beam coordinates with origin at tip 

Dimensionless coordinates of beam tip in X,Y system 

Angle between tangent to beam neutral axis and 
horizontal (slope of neutral axis) 

Angle between tip load and vertical 

Angle between horizontal and line between root and tip 

BEAM EQUATION 


The beam is considered horizontal, clamped at the right 
end and loaded at the left end as shown in figure 1. The 
coordinate system for the beam is shown in figure 2. The 
origin of the coordinates X,Y is at the clamped end or 
root, whereas the origin of coordinates XI, Y1 is at the 
tip. The beam is considered i nextensi onal , and transverse 
shear deformations are neglected. Equilibrium equations are 
written with reference to the element shown in figure 3. 
Horizontal equilibrium gives 

(?+ dF)cos(4>4-ci^) — i|> -V-Qs\\r> vp — ($+ <AQ) c4\p) O 



By performing algebraic manipulations and neglecting terms 
containing products of differentials, there results 

4> — CQ — comstanV - P sWi 

Vertical equilibrium gives 

(P-V ap') tc* i{>) + CQ-V do) C 05 (i\» V a 4>) *” Ps'in 4> - Q c*> *, = o 

Similar manipulations lead to 

P-Vlvn^J -V Qcasxty — P cos 

Equilibrium of moments gives 

tA- N\- drv\ - QdlS =.0 

or 

4K 

The load P along the neutral axis of the beam can be 
eliminated from equations (2) and (4) to give 

Qc p -V 4©S T? Co s 


Engineering beam theory gives 

Substituting equations (6) and (7) into equation <5) results 
in the nonlinear beam differential equation 

- — ~~0bun if CoS f<U> f, ^>) 

d€> ^ " 


Equation (8) can be made dimensionless by simply defining 



so that equation (8) becomes 
or 

f* jk - - a. 'V 4- t ^ 40 

where the quantity A=Fl?/EI is termed, in this paper, the 
1 oad-st i f f ness parameter. 


NUMERICAL INTEGRATION OF BEAM EQUATION 

Equation (9) can be integrated numerically in a direct 
fashion. The second-order differential equation is replaced 
by two first-order differential equations 



du 

dS 


=• - A (s-’w'i £ w V-f tcs £ cos\/) 


dV _ 
ds 


o 


Equations (10) are integrated 
Euler method (ref. 3). Let 


numerically with the modified 


-f V) = - A (*»* ^ sih V -v C 05 . £^sV) 

If the symbol h denotes the increment in the length 
parameter s , and the subscripts 0 and 1 denote values 
at the beginning and end of an increment, respectively, then 


UjC\) =. U 0 -V h £ ($ } Vo) 

V,(l) = Vg + h u 0 

U, n> = Uo + 1 W Cf ,V„) + +<■$, v, (0)1 

V, «> = M 0 4 ^LUe * 0,<>)] 

u,Bl = u 0 + £[f L^Yo) + f Of, V, <z))] 

V,W = V D + h-[U 0 + U,«)1 


00 


V (V2) 


U,CN> = U 0 4^Lf«V 0 )+^4, V.CK-0X1 

n/.WV = v 0 -v|iLu o -vo t (M-0] 


The quantity h is assumed constant over the length of the 
beam. The accuracy of the calculations depends on h and is 
improved by reducing the sice of h . The numbers 
(1,2,...,N> in parentheses in equations (12) denote various 
approximations to the values of U and V at the end of the 
increment. Continue the iteration until successive values of 
U* or are equal. Then replace U 0 by L) t and V 0 by Oj 

and advance to the next increment. This integration process 
is started at the tip of the beam by assuming a value of the 
slope of the beam V 0 and setting the curvature * U o =0 . The 
integration continues toward the root until the slope becomes 
zero. The whole process requires iteration because the value 
of the load-stiffness parameter A must be adjusted to make 
the slope zero at the end of the final integration increment. 
(See ref. 3 for a related approach to the cantilever beam 
with uniform distributed load.) 

A BASIC computer program which integrates equations (10) 
using the process laid out in equations (12) is presented in 
the Appendix. With practice, given starting values at the 
tip for beam slope? V 0 and beam curvature U o =0 , the user of 
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this program can converge in just a -few iterations to values 
of A which lead to zero slope at the root of a beam of unit 
1 ength . 


RESULTS 


Results are presented for four load cases: (1) tip load 

vertical (vertical load), (2) tip load normal to a straight 
line connecting root and tip (slewing load), (3) tip load 
normal to the beam (follower load), and (4) tip load 
horizontal (column load). All results presented in the paper 
are calculated with the increment h=0.01 . 

For cases (1) vertical load, (3) follower load, and (4) 
column load, the quantity t, is simply prescribed initially, 
and the calculation follows directly. For case (2) slewing 
load, the quantity £ must be equal to & . Since & is not 

known a priori, an additional tr i al -and-error operation is 
required to make £ equal to 0 . 

Results are presented in Tables 1 and 2 and in figures 
4-9. Values of the load-stiffness parameter A and the 
angle between the horizontal and the line between root and 
tip of the beam 0 are listed in Table 1 for various values 
of tip slope. Also in Table 1 are comparisons with results 
of closed form solutions given in references i and 4. The 
curves in figures 4-7 represent dimensionless deflection 
shapes of beam neutral axes. Tip load vectors are shown in 
figures 4-7, and these vector magnitudes are indeed 
proportional to the tip load magnitudes. 

For the case of vertical tip load, the load-stiffness 
parameter is plotted as a function of dimensionless vertical 
tip deflection in figure 8. The contrast between linear and 
nonlinear solutions for this case is striking. For the case 
of the slewing tip load, the load-stiffness parameter is 
plotted in figure 9 as a function of Q , the angle between 
the horizontal and the line between root and tip. For & up 
to one radian the beam behavior is nearly linear in this 
slewing load case. 


CONCLUDING REMARKS 


Equations are presented for the numerical solution of the 
nonlinear, large deflection beam equation for cantilever 
beams under tip loads of various orientations. Results are 
presented for selected cases, and comparisons are made with 
closed form solutions from the literature. The numerical 
solution, implemented in a BASIC-1 anguage computer program, 
is accurate. The computer program can operate on 
microcomputers, and solutions can be obtained rapidly and 
cheaply. 
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APPENDIX 


COMPUTER PROGRAM FOR LARGE DEFLECTION OF CANTILEVER BEAM 
UNDER ARBITRARILY DIRECTED TIP LOAD 


SYMBOL IN TEXT 

VARIABLE. IN COMPUTER PROGRAM 

A 

A 

El 

El 

F 

F 

h 

H 

L 

L 

s 

SBAR 

>° 

r. 

D° 

UO, VO 

u, 

Ui , VI 

XI, Y1 

XI, Y1 

XL, YL 

XL, YL 


PS I 

e 

THETA 

i 

ZETA 



ulist 


10 REM INTEGRATION OF BEAM EQUATIONS 021985 

20 REM LARGE DEFLECTION OF CANTILEVER BEAM UNDER ARBITRARILY DIR 
ECTED TIP LOAD 

30 REM BASIC DIFFERENTIAL EQUATION: D"'2PSI /DSBAR / ‘'2=-A* (SIN (ZETA) * 
SIN (F'SI ) +COS <2 ETA) *COS (PS I ) ) 

40 REM PSI=SLOPE OF BEAM: SBAR=S/L: S=DISTANCE ALONG BEAM: L=LENGTH 
OF BEAM: A= (F*L-''2> / (E*I > : ZETA=ANGLE OF TIP LOAD WITH VERTICAL 
50 REM WRITE SECOND ORDER DIFFERENTIAL EQUATION AS TWO FIRST ORD 
ER DIFFERENTIAL EQUATIONS 

60 REM DU/DSBAR=— A* (SIN ( ZETA) *SIN (V) +COS (ZETA) *COS ( V) ) 

70 REM U=DV/DSBAR 

80 REM U=DPSI /DSBAR: V=PSI 

85 INPUT "TODAYS DATE= " ; D 

90 INPUT "SLOPE AT TIP IN RADI ANS= " j VO 

100 INPUT "CURVATURE AT TIP= " j UO 

110 INPUT "ANGLE OF TIP LOAD WITH VERTICAL IN RADI ANS= " j ZETA 
120 INPUT "INCREMENT IN LENGTH PARAMETER SBAR= ";H 
130 INPUT "LOAD-STIFFNESS COEFFICIENT= ";A 

135 INPUT "DO YOU WANT TO PRINT X AND Y? 1 IF YES, 2 IF NO: ";Q1 

140 PRINT CHR$ <4>;"PR#1": PRINT CHR$ (9) j "80N" 

145 PRINT "RESULTS FROM BEAM I NTEGR AT I ONT I PL0AD02 1 985 CALCULATED 0 ; 
n " ; D I 

150 PRINT "SLOPE AT TIP IN RADI ANS= " ; VO 
160 PRINT "CURVATURE AT TIP= " ; UO 

170 PRINT "ANGLE OF TIP LOAD WITH VERTICAL IN RADI ANS= " ; ZETA 
180 PRINT "LOAD-STIFFNESS COEFFICIENT^ ";A 

190 PRINT "INCREMENT IN LENGTH PARAMETER SBAR= ";H 

191 IF Q1 < >1 THEN 200 

192 PRINT : PRINT : PRINT 
194 POKE 36,5: PRINT "XI"; 

196 POKE 36,25: PRINT "Yl" 

200 PRINT CHR$ (4) 5 "PR#0" 

210 REM Ul < 1 ) =UO+H* (-A) * (SIN < ZETA) *SIN (VO) +COS ( ZETA) *COS (VO) ) 

220 REM VI (1)=V0+H*U0 

230 REM FOR N>1 , Ul <N) =UO+H* (-A) * (SIN ( ZETA) *SIN (VO) +COS ( ZETA) *COS 

(VO) +SIN (ZETA) *SIN (VI (N-l) ) +COS ( ZETA) *COS (VI (N-l) > ) /2 

240 REM FOR N>1, VI (N) =VO+H* (UO+Ul (N-l > ) /2 

250 REM ITERATE UNTIL SUCCESSIVE VALUES OF VI ARE EQUAL 

260 REM THEN REPLACE UO AND VO BY Ul AND VI, RESPECTIVELY 

270 REM CONTINUE UNTIL V1<0 

280 REM A1=U1 (1) :A2=U1 (2) :B1=V1 (1) :B2=V1 (2) 


290 XI 

1! 

t a 

< 

= 0 : c = 

0:J = 

1 







295 DIM XI (10) , Yl (10) 









300 A1 

= UO + 

H * ( - 

A) * ( 

SIN 

(ZETA) 

* SIN 

(VO) 

4- 

COS 

(ZETA) 

* COS 

(VO) ) 










310 B 1 

= VO + 

H * UO 









320 A2 

= UO + 

H * ( - 

A) * ( 

SIN 

(ZETA) 

* SIN 

(VO) 

+ 

cos 

(ZETA) 

* COS 

( VO ) + 

BIN (ZETA) * 

SIN 

(Bl) + 

COS (ZETA) 

* 

cos 

(Bl) ) 


330 B2 = VO +H* (U0+A1) /2 
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335 POKE 36,0: PRINT A2; 

336 POKE 36,20: PRINT B2 

340 IF B2 = B1 THEN 1000 
350 A1 = A2 

360 B1 = B2 

370 GOTO 320 

1000 POKE 36,0: PRINT A2; 

1010 POKE 36,20: PRINT B2 
1015 PRINT 

1020 IF B2 < 0 THEN 1090 
1030 XI = XI + H * COS (B2) 

1040 Y1 = Y1 + H * SIN <B2) 

1044 C = C + 1 

1045 IF C = (.1 / H) * J THEN 3000 

1050 UO = A2 

1060 VO = B2 

1080 GOTO 300 

1090 F< = (VO * H) / (VO + ABS (B2> ) 

1095 XI = XI + R * COS (VO) 

1100 Yi = Y1 + R * SIN (VO) 

2000 PRINT CHR$ (4) j "PR#1 " : PRINT CHR$ (9) ; "SON" 
2010 PRINT . 

2020 POKE 36,0: PRINT "XL= ";X1; 

2030 POKE 36,20: PRINT ” YL= ";Y1 
2040 THETA = ATN (Yl / XI) 

2050 POKE 36,0: PRINT "THETA= THETA; 

2060 L = H * C + R 

2070 POKE 36,20: PRINT "DISTANCE TO ZERO SLOPE= " ; L 

2085 IF (31 < > 1 THEN 2565 

2090 GOTO 2450 

2450 PRINT : PRINT : PRINT 

2460 POKE 36,5: PRINT "X"; 

2470 POKE 36,25: PRINT "Y" 

2500 PRINT 

2502 POKE 36,0: PRINT XI 5 

2504 POKE 36,20: PRINT Yl 

2505 FOR J = 1 TO 10 
2510 X = XI (10) - XI (J) 

2520 Y = Yl (10) - Yl ( J ) 

2530 PRINT CHR* (4)J"PR#1 H : PRINT CHR* (9) ; "SON" 
2540 POKE 36,0: PRINT X; 

2550 POKE 36,20: PRINT Y 
2560 NEXT J 

2565 PRINT CHR$ (4) ; "PR#0" 

2570 END 

3000 IF Q1 < > 1 THEN 1050 

3010 PRINT CHR* (4);"PR#1": PRINT CHR* (9) j "BON" 
3030 POKE 36,0: PRINT XI; 

3040 POKE 36,20: PRINT Yl 
3050 XI (J) = XI 
3060 Yl ( J ) = Yl 
3070 J — J + 1 

3090 PRINT CHR$ (4) j "PR#0" 

3100 GOTO 1050 
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TYPICAL PRINTOUT 


RESULTS FROM BEAM I NTEGRAT I ONT I PLO AD02 1 9S5 CALCULATED ON 52085 
SLOPE AT TIP IN RADI ANS= 1.3962634 
CURVATURE AT TIP= 0 

ANGLE OF TIP LOAD WITH VERTICAL IN RADI ANS= 1.57079633 
LOAD-STIFFNESS COEFFICIENT= 3.19267 
INCREMENT IN LENGTH PARAMETER SBAR= .01 


XI 

.0179601298 
. 039 1 506505 
. 0665560747 
. 102929863 
. 150564535 
.21096228 
.234447718 
.369330612 
.464276777 
. 563530959 


Y1 

. 0983726294 
. 196091612 
. 2922358 
. 38533243 
.473167718 
.552725781 
. 620330257 
. 672040657 
. 704279598 
.714550956 


XL— .563531231 YL.= .714550956 

THETA= .903014 DISTANCE TO ZERO SLOPE= 1.00000027 
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TYPICAL PRINTOUT (Continued) 



X 

.563531231 
. 545570829 
. 524380308 
.496974884 
. 460601095 
.412966423 
. 352568679 
.279083241 
. 193700347 
.0992541814 
0 


Y 

.714550956 
.616178327 
.518459344 
.422315156 
.329218526 
. 241 3832.^8 
. 161825176 
. 0942206995 
.042510299 
.0102713581 
0 
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TABLE 1. RESULTS ON LARGE DEFLECTION OF C .ANTI LEVER BEAMS 
UNDER VARIOUS TIP LOADS 


Tip Slope 


Theta F/ (Buckling Load) 


Radians 

(This paper) 

(Ref. 1) 

Radians 

(This paper) 



0. 17453293 
0.34906585 
0.52359878 
0.69813170 
0.87266463 
1.04719755 
1 . 22173048 
1 . 39626340 





VERTICAL TIP LOAD 


0 . 353 

0. 731 

1. 163 
1 . 692 

2 . 7192 

3 . 405 
5.081 
8.678 



3.405 

8.678 


0. 116 

0. 143 

U ■ 2o*2 

0.296 

0.349 

0.471 

0.469 

0.686 

0.592 

0.970 

0.722 

1 . 380 

0.861 

2.059 

1 . 024 

3.517 


SLEWING TIP LOAD 


0. 17453293 

0.349 

0. 34906585 

0. 701 

0.52359878 

1 . 057 

0. 69813170 

1 . 420 

0.87266463 

1.791 

1.04719755 

2. 174 

1 . 22173048 

2.570 

1 . 39626340 

n OQ r ? 



FOLLOWER TIP LOAD 


0. 17453293 

0.349 

0.34906585 

0. 701 

0.52359878 

1 . 057 

0. 69813170 

1 . 420 

0. 87266463 

1 . 791 

i. 04719755 

2. 175 

1.22173048 

2.575 

1 . 39626340 

2.994 



O. 115 
0.231 
0. 346 
0. 461 
0.575 
0.689 
O. 801 
O. 913 


0. 17453293 

2.477 

0. 34906585 

2.505 

0. 52359878 

2.554 

0.69813170 

2.624 

0.87266463 

2.719 

1 . 04719755 

2.842 

1.22173048 

2.997 

1 . 39626340 

3. 193 


HORIZONTAL TIP LOAD 


O. 1 10 
0.221 
0.332 
0. 443 
O. 556 
0 .670 
0.785 
0 . 903 



0. 142 
0.284 
0.428 
0.575 
0.726 
0.882 
1.044 
1.214 


1 . 004 
1.015 
1 . 035 
1 . 064 
1 . 102 
1. 152 
1.215 
1.294 




0.116 

0. 142 

0 .231 

0.284 

0.347 

0.428 

0.463 

0.575 

0. 580 

0.726 

0.698 

0.881 

0.816 

1 . 042 

0.935 

1 . 209 





1.015 


1 . 063 
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TABLE 2. TIP COORDINATES FROM LARGE DEFLECTIONS OF 
CANTILEVER BEAMS UNDER VARIOUS HP LOADS 


Up Slope 


Degrees 


Dimensionless Tip Coordinates 


(This paper) 


(Ref. 4) 


VERTICAL TIP LOAD 


10 

0.992 

. 0. 115 

20 

0. 968 

0.228 

30 

0.928 

0.338 

40 

0.873 

0.442 

50 

0.802 

0.540 

60 

0.716 

0.630 

70 

0.611 

0.712 

80 

0.481 

0. 790 

SLEWING T] 

10 

0.992 

0. 115 

20 

0.968 

0.228 

30 

0.929 

0. 336 

40 

0.875 

0. 437 

50 

0.808 

0.529 

60 

0.729 

0.611 

70 

0.640 

0.680 

80 

0. 543 

0.735 





0. 992 
O. 963 
0.929 
0. 876 
0.810 
O. 734 
0.647 
0. 557 


FOLLOWER TIP LOAD 


0. 115 
0.228 
0 . 335 
0.435 
0. 525 
0. 604 
0.672 
0.721 



HORIZONTAL TIP LOAD 



0.992 

0. 110 

0.970 

0.218 

0.933 

0.321 

0.882 

0.419 

0.819 

0.509 

0. 744 

0.589 

0 . 658 

0.658 

0.564 

0.715 


0.970 

0.220 

0. 881 

0.422 

0.741 
0. 560 

0.593 

0.719 
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Figure 5. Deflection curves for cantilever beam under slewing tip load. 
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Figure 9. Load-stiffness parameter as a function of angle 0. (between 

horizontal and line from root to tip) for cantilever beam under 
slewing tip load. 
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